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ABSTRACT 

We report a study of the spectra and structure of the accretion disk of the dwarf nova V2051 
Ophiuchi while the star was in an unusual faint brightness state during 1996. The differences 
between the UV resonant lines and continuum disk surface brightness distributions indicate a 
vertically extended disk with the emission from these lines arising from the upper atmospheric 
layers. Distinct emission along the stream trajectory suggests the occurrence of gas stream 
overflow. Spatially resolved spectra show that the lines are in emission at all disk radii. The 
Balmer decrement becomes shallower with increasing radius. The FWHMs of the emission lines 
show differences with respect to that expected for a gas in Keplerian rotation and the line 
intensities drop with a radial dependency of / oc The uneclipsed light contributes about 

5-10 per cent of the total flux, and its spectrum is dominated by a Balmer jump and strong lines 
in emission. Broad absorption bands, possibly due to Fell, are present in the spectra of the disk 
side farther away from the secondary star, suggesting it arises from absorption by a extended 
gas region above the disk; the differences between the spectra of the hemisphere farther from 
and nearer to the secondary star are interpreted in terms of chromospheric emission from a disc 
with a non-negligible opening angle (limb brightening effect). Stellar atmosphere model fits to 
the extracted white dwarf spectrum lead to a temperature Two — SSOO^^annif and a distance 
of d = 67 if the inner disk is opaque, or d 



'1900 ' 

92 isspc if the inner disk is optically thin. 
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1. Introduction 

V2051 Ophiuchi was discovered by Sanduleak 
(1972). It is a dwarf nova, a compact binary in 
which mass is fed to a white dwarf (the primary 
star) by a Roche lobe filling late-type star (the sec- 
ondary) via an accretion disk. V2051 Oph shows 
recurrent but sparse outbursts in which the disk 
expands and brightens by 2-3 magnitudes during 
1-3 days (B ~ 13 mag at maximum, Bateson 1980; 
Warner & Cropper 1983; Warner & O'Donoghue 
1987) as a consequence of a sudden increase in 
mass accretion through the disk. In the quiescent, 
low-mass accretion state, the white dwarf and the 
bright spot (formed by the impact of the infalling 
gas stream with the outer edge of the disk) dom- 
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inate the light from the binary at optical and ul- 
traviolet wavelengths, and the optical spectrum 
shows double-peaked HI and Hell emission lines 
which exhibit the classical rotational disturbance 
effect during eclipse (Cook & Brunt 1983; Watts 
et al. 1986). 

The binary is seen at a high inclination angle 
{i = 83°), which leads to deep eclipses {AB ~ 
2.5 mag) in the light curve every 1.5 hours when 
the white dwarf, accretion disk and bright spot 
are occulted by the mass-donor star. This al- 
lows the emission from the different light sources 
to be distinguished and spatially resolved stud- 
ies to be performed, making V2051 Oph an excel- 
lent laboratory for the study of accretion physics. 
In particular, it is well suited for the application 
of indirect imaging techniques to resolve the disk 
emission both in position (eclipse mapping, Horne 
1985) and velocity (Dopplcr tomography, Marsh 
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& Home 1988). 

V2051 Oph displays large amplitude flickering 
activity (up to 30 per cent of the total light in 
the optical), which is responsible for a variety 
of eclipse morphologies (e.g., Warner & Cropper 
1983) and usually hampers the measurement of 
white dwarf and bright spot eclipse phases. Bap- 
tista ct al. (1998) caught the star in an unusual 
faint brightness state (B ~ 16.2 mag) in which 
the mass accretion rate and flickering activity were 
significantly reduced, allowing a clean view of the 
white dwarf at disk center. They derived a mass 
ratio of <7 = 0.19 ± 0.03 and an inclination of 
i = 83.°3± l.°4. 

Vrielmann et al. (2002) applied eclipse mapping 
techniques to multicolor data of V2051 Oph to find 
that the disk is brighter in the front side (the hemi- 
sphere nearer to the secondary) and interpreted 
this behavior as caused by enhanced emission from 
the bright spot. They also estimated a distance of 
d = 146 ± 20 pc for the system. 

Kiyota & Kato (1998) observed V2051 Oph in 
supcroutburst and detected supcrhumps in the 
light curve of the object, confirming its classifi- 
cation as an SU UMa star - a sub-class of dwarf 
novae with occasional supcrout bursts ^ 0.7 mag 
brighter and ~ 5 times longer than an ordinary 
outburst. However, due to the complex behavior 
observed in V2051 Oph, Warner (1996) suggested 
a possible alternative classification of it as a po- 
laroid, an intermediate polar with a synchronized 
primary. 

Here we report the results of an optical and UV 
spectral mapping experiment of V2051 Oph in its 
faint brightness state of 1996 (e.g. Baptista et al. 
1998). The details of the data analysis are given 
in Section 2. In Section 3 we investigate the disk 
spatial structure in the lines and continuum, and 
we present spatially resolved spectra of the disk, 
gas stream and uneclipsed light. The results are 
discussed in Section 4 and summarized in Section 
5. 

2. Data analysis 
2.1. Observations 

The Faint Object Spectrograph (FOS) onboard 
the Hubble Space Telescope (HST) was used to 
secure high-speed spectroscopy covering two con- 



secutive eclipses of V2051 Oph on 1996 January 
29 while the star was in an unusual low brightness 
state {B ~ 16.2 mag). The observations are sum- 
marized in Table 1. The first nm covers a short 
phase range around the egress of an eclipse in the 
optical (G400H grating), and the second run cov- 
ers the next eclipse in the UV (G160L grating). 
The reader is referred to Baptista et al. (1998) for 
a detailed description of the observations and data 
reduction procedures. 

The average out-of-eclipse UV spectrum is 
shown in Fig. 1 (top). It displays emission lines of 
Lya A1216 (mostly geo coronal), N V A1240, 1243, 
Sin A1300, CII A1336, Si IV A1394, 1403, CIV 
A1549,1551, SilH A1892, CIII A2297, as well as 
broad absorption bands reminiscent of those seen 
in OY Car (Horne et al. 1994), which were in- 
terpreted as a blend of Fell lines. There is no 
evidence of Hell A1640 emission. The position 
of the second-order Lya emission is illustrated in 
Figure 1 (top). The bottom panel of Figure 1 
shows the average out-of-eclipse optical spectrum. 
It displays the Balmer continuum in emission, as 
well as He I, Hell and strong Balmer emission 
lines. The optical emission lines show a clear 
double-peaked profile. The blue peak of the lines 
is stronger than the red peak because part of the 
disk side contributing redshifted line emission is 
still occulted by the secondary star at the phase 
range used to compute the average spectrum. 

2.2. Light-curve construction 

The UV spectra were divided into a set of 34 
narrow passbands, with 22 continuum passbands 
19-60 A wide and 12 passbands for 10 lines (Fig. 1). 
The blue end of the spectra (shortward of the Lya 
line) was not included in this analysis because of 
the very low count rate [and, consequently, low 
signal-to-noise ratio (S/N)] of the data in this 
spectral region. The emission lines were sampled 
in a single velocity bin 5000 or 6000fcms~"'^ wide 
(centered at the rest wavelength, v = Okms~^), 
except for the C IV line, which was separated in a 
bin centered at the rest wavelength and velocity- 
resolved bins for the blue and red wings of the 
line (all bins with Av = 5000fcms~^)^ The opti- 

Thc choice of the number (or width) of the passbands was 
constrained by the signal-to- noise (S/N) ratio of the result- 
ing light curves. Narrower passbands lead to light curves 
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Fig. 1.— Average out-of-eclipse UV (G160L; top ) and optical (G400H; bottom) HST spectra of V2051 Oph. 
The phase ranges arc +0.10 to +0.20 and +0.075 to +0.080 cycles, respectively, for the UV and optical 
spectra. The observations correspond to the eclipse cycles 109,988 (optical) and 109,989 (UV) according to 
the linear ephemeris of Baptista et al. (2003). Emission (vertical dotted lines) and absorption (horizontal 
lines) features are labeled. Horizontal tick marks indicate the passbands used to extract light curves (34 in 
the UV and 68 in the optical). 



Table 1: Journal of the Observations. 



Run Start Eclipse Number of Spectral Range Phase Range 

(UT) Cycle Exposures A (cycles) 

"HI 19!03 109,988 US 3226 - 4781 +0.01, +0.08 

H2 20:24 109,989 693 1150-2507 -0.09, +0.34 
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cal spectra were divided into a set of 68 passbands, 
with 32 continuum passbands 15-42 A wide and 36 
passbands for 11 lines. The emission hnes were re- 
solved in velocity bins 1000 km wide, with one 
bin centered at the rest wavelength. The systemic 
velocity of V2051 Oph (7 ~ 40fcms-i; Watts 
et al. 1986) is both rather uncertain and much 
smaller than the width of the passbands and was 
neglected. For those passbands including emission 
lines the light curve comprises the total flux at the 
corresponding bin, with no subtraction of a possi- 
ble continuum contribution. 

Light curves were extracted for each passband 
by computing the average flux on the correspond- 
ing wavelength range and phase-folding the results 
according to the linear ephemeris of Baptista et al. 
(2003), 

Tmid(BJDD) = 2 443 245.97752+0.062 427 8634 E, 

(1) 

where Tmid gives the inferior conjunction of the 
white dwarf. 

Light curves at two selected continuum pass- 
bands and for the C IV and H7 lines can be seen in 
Figs. 2 and 3, respectively. The incomplete eclipse 
phase coverage of the HI run is clear (Fig. 3). 

2.3. Eclipse mapping 

Maximum-entropy eclipse mapping techniques 
(Horne 1985, Baptista & Steiner 1993) were used 
to solve for a map of the disk brightness distri- 
bution and for the flux of an additional uneclipsed 
component in each band. The reader is referred to 
Baptista (2001) for a recent review on the eclipse 
mapping method. 

As our eclipse map we adopted a flat grid of 
51 X 51 pixels centered on the primary star with 
side 2i?ii, where Rli is the distance from the 
disk center to the inner Lagrangian point. The 
eclipse geometry is defined by the inclination i and 
the mass ratio q. The mass ratio q defines the 
shape and the relative size of the Roche lobes. The 
inclination i determines the shape and extension of 
the shadow of the secondary star as projected onto 
the orbital plane. In this paper we adopted the 
values derived by Baptista et al. (1998) , q = 0.19 
and i = 83°, which correspond to an eclipse phase 

too noisy to produce reliable eclipse maps. 



width of Acf) ~ 0.0662 cycle. This combination of 
parameters ensures that the white dwarf is at the 
center of the map. 

For the reconstructions we adopted the default 
of limited azimuthal smearing of Rutten et al. 
(1992), which is better suited for recovering asym- 
metric structures than the original default of full 
azimuthal smearing (e.g. Baptista et al. 1996). 
We used a radial blur width Ar = 0.0157i?Li and 
an azimuthal blur width = 30°. 

The statistical imcertainties in the eclipse maps 
were estimated with a Monte Carlo procedure 
(e.g., Rutten et al. 1992). For a given narrow band 
light curve, a set of 20 artificial light curves was 
generated in which the data points were indepen- 
dently and randomly varied according to a Gaus- 
sian distribution with standard deviation equal to 
the uncertainty at that point. The light curves 
were fitted with the eclipse-mapping code to pro- 
duce a set of randomized eclipse maps. These 
were combined to produce an average map and 
a map of the residuals with respect to the av- 
erage, which yields the statistical uncertainty at 
each pixel. The uncertainties obtained with this 
procedure were used to estimate the errors in the 
derived radial brightness distributions, as well as 
in the spatially resolved spectra. 

The Appendix presents the results of eclipse 
mapping simulations addressing the reliability of 
disk surface brightness reconstructions from light 
curves of incomplete phase coverage, such as those 
of run HI. 

Light curves and respective eclipse maps at se- 
lected passbands for the continuum and for the 
strongest emission lines are shown in Figs. 2 and 
3, respectively. Contour curves in these Figures 
enclose the regions of the eclipse maps at the 3 a 
and 5 a levels of statistical significance. 

3. Results 

3.1. Disk structure 

Maps of the disk surface brightness distribu- 
tions, calculated by the maximum entropy eclipse 
mapping method, allow us to obtain spatially re- 
solved information about the disk emission. In 
this section we discuss the structures in the eclipse 
maps of the strongest lines in the spectrum, as well 
as in selected continuum passbands. 
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Fig. 2. — Data and model light curves and eclipse maps for the C IV line center (Aw — 5000 km s~^) and for 
two selected continuum passbands in the UV. Left: Data (dots with error bars) and model (solid lines) light 
curves. Dotted lines mark the contact phases of the white dwarf eclipse. Right: Corresponding eclipse maps 
in a logarithmic gray scale. Brighter regions are shown in black, fainter regions in white. Contour curves 
enclose the eclipse map regions at the 3 a and 5 a levels of statistical significance. A cross marks the center 
of the disk; dotted lines show the Roche lobe, the gas stream trajectory, and a disk of radius 0.56 Rli', the 
secondary is to the right of each map and the stars rotate counter-clockwise. The vertical bars indicate the 
logarithmic intensity level of the gray scale in each map. 
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Fig. 2 shows the hght curves and respective 
echpse maps for the CIV hne center and for two 
selected UV continuum passbands. The S/N ra- 
tio of the hght curves decreases towards the blue 
end of the spectrum. Accordingly, the statisti- 
cal significance of the eclipse maps is lower in this 
spectral region. Thus, the brightest parts of the 
1444 A continuum map are only significant at the 
3(7 level, while the eclipse maps in the red side 
of the spectrum are typically significant at the 5 cr 
confidence level. 

The CIV line center light curve shows an in- 
verse "Gaussian" -shaped eclipse with a minor 
asymmetry at egress and no sharp breaks in 
slope. Accordingly, the eclipse map shows a broad, 
smooth, and fairly symmetrical brightness distri- 
bution centered on the white dwarf, with no clear 
sign of the bright spot. 

In contrast to what is observed in the CIV 

line center, the continuum light curves show sharp 
breaks in slope at the ingress and egress phases 
of the white dwarf plus a conspicuous asymmet- 
ric egress shoulder, indicating that there is signif- 
icant additional emission from the disk side con- 
taining the bright spot. The corresponding eclipse 
maps display pronounced emission from the com- 
pact white dwarf at disk center plus an asymmet- 
ric structure in the disk side that is moving away 
from the secondary star (the upper hemisphere of 
the eclipse maps in Fig. 2) that can be associated 
with the gas stream emission. 

The presence of bright compact structures in 
the continuum maps and their absence in the C IV 
line map, in addition to the observed spectrum of 
the uneclipsed component (see Section 4.4), are 
indications of the vertical extension and large op- 
tical depth of the gas from which the CIV line 
originates. Our interpretation is that the C IV line 
emission arises in a large gas region surrounding 
the disk. This region needs to be vertically ex- 
tended and optically thick at this wavelength in 
order to veil the strong emission from the under- 
lying white dwarf and bright spot. This is in line 
with expectations for an emission line produced 
by resonant scattering in a disk chromosphere. A 
similar result was found for the nova-like variable 
UX UMa (Baptista et al. 1995). 

The eclipse map at 2251 A is representative of 
most of the continuum maps. Aside from the 
asymmetry caused by the bright spot and gas 



stream emission, it shows that another brightness 
asymmetry, namely, the disk hemisphere closest 
to the secondary star (the "front" disk side) is 
systematically brighter than the disk side farther 
away from the secondary star (hereafter called the 
"back" disk side). A similar effect was observed 
by Vrielmann et al. (2002), who interpreted it as 
a consequence of azimuthally-smeared bright spot 
emission. Our results will lead to a different inter- 
pretation (see Section 4.1). 

Fig. 3 shows light curves and eclipse maps for 
the if 7 line center {v = Okms~^) and the blue 
{v = — 1000km s~^) and red {v = -1-1000 fcm s^^) 
line wing passbands. The velocity-resolved Hj 
line light curves show the expected behavior for 
the eclipse of gas rotating in the prograde sense 
(rotational disturbance), with the blue wing of the 
line reappearing from eclipse earlier than the red 
wing. The phase of maximum derivative in the 
three light curves coincides with the white dwarf 
egress phase, indicating that the maxima of the 
distributions are at the disk center. The light 
curves show deep eclipses with smooth and long 
lasting egresses. As a consequence, the resulting 
maps show fairly broad brightness distributions 
centered at the white dwarf position with an asym- 
metric bright source in the quadrant containing 
the bright spot. Simulations show that, in spite of 
the incomplete phase coverage of the optical light 
curves, it is possible to derive a fairly good recon- 
struction of the position and intensity of the white 
dwarf and bright spot in the corresponding eclipse 
maps (see the Appendix). 

Fig. 3 (third and fourth columns) shows the net 
line emission and the net line absorption in each 
case. Net line emission/absorption maps are ob- 
tained by combining continuum eclipse maps on 
the short- and long-wavelength sides of the tar- 
get emission line and by subtracting the derived 
average continuum map from each of the velocity- 
resolved line maps. Positive intensities in the re- 
sulting maps (third column, black and gray) sig- 
nal the regions where the line appears in emission, 
while negative intensities (fourth column, black 
and gray) trace the regions where the line is in 
absorption. The net emission/absorption maps 
reveal that the Balmer lines are in emission over 
most of the accretion disk, but appear in absorp- 
tion at the disk center and also in the region of 
the bright spot at the disk rim. 
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0.025 0.05 0.075 0.1 total emission net line emission net line absorption 
orbital phase (line + continuum) 



Fig. 3. — Left: Data (dots with error bars) and corresponding model (solid lines) light curves for the line 
center passband {v = Ofcms"^, Aw = 1000 fcms~^) and for its blue and red line wings (Au = 1000 fcms~^). 
Second from left: Corresponding eclipse maps (total emission). Third from left: Eclipse maps of the net line 
emission. Right: Eclipse maps of the net line absorption in each case. The gray scale is the same in each 
column. The notation is similar to that of Fig. 2. 



7 



3.2. Spatially resolved spectra 

Our set of monochromatic eclipse maps allows 
us to separate the spectra emitted by different 
parts of the accretion disk. Motivated by the dis- 
tinct emission observed in the gas stream region 
and by the systematic difference in the emission 
from the opposite disk hemispheres, we sliced the 
disk into three distinct regions: "front", "back" 
and "gas stream". Azimuths are measured from 
the line joining both stars and increase counter- 
clockwise. We define "front" as the disk section 
between 270° and 90°, and "back" as the region 
between azimuths 90° and 270°. The region de- 
fined as "gas stream" is depicted in Fig. 4 (gray). 
In order to separate the disk spectra at different 
distances from the disk center, we further divided 
each of the disk regions in concentric annular sec- 
tions of width AR = 0.1i?Li . 

In order to minimize the possible contribu- 
tions from the bright spot and gas stream to the 
disk "front" and "back" spectra we calculated the 
symmetric disk-emission component in each hemi- 
sphere. The symmetric component is obtained by 
slicing the disk into a set of radial bins and by fit- 
ting a smooth spline function to the resulting set 
of medians of the lower quartile of the intensities in 
each bin. The spline fitted intensity in each annu- 
lar section is taken as the symmetric component. 
This procedure essentially preserves the baseline 
of the radial profile, removing all azimuthal struc- 
ture. The statistical uncertainties affecting the fit- 
ted intensities are estimated with the Monte Carlo 
procedure described in Section 2.2. 

The spatially resolved disk (front and back) 
spectra are shown in Fig. 5. The inner disk shows 
a flat continuum which becomes progressively 
fainter and redder with increasing radius, indicat- 
ing the existence of a radial temperature gradient. 
The lines and the Balmer jump are in emission 
at all disk radii. The Balmer decrement becomes 
flatter and the lines are more prominent and nar- 
rower with increasing radii. The UV spectrum of 
the "front" side is perceptibly different from that 
of the "back" side. The comparison of the UV 
"front" and "back" spectra shows that the latter 
contains broad absorption bands, possibly due to 
Fell, that become more conspicuous in the outer 
parts of the disk - suggesting that it arises from 
absorption along the line of sight by a vertically 
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Fig. 4. — Schematic diagram showing the re- 
gions defined as "front" , "back" and "gas stream" . 
Dashed lines mark the annular regions of width 
0.1 used to extract spatially resolved spectra. 
Dotted lines show the projection of the primary 
Roche lobe onto the orbital plane and the gas 
stream trajectory. Azimuths are measured with 
respect to the line joining both stars and increase 
counter-clockwise. Four reference azimuths are la- 
beled. 
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Fig. 5. — UV (left) and optical (right) spatially resolved disk spectra for the azimuthal sections defined as 
front (270° — 90°) (gray) and back (90° — 270°; black lines) for three selected annular regions (labeled in 
units of Rli)- Vertical dotted lines mark the major transitions of the spectra. 



9 




1500 2000 2500 3500 4000 4500 

wave length (A) 



Fig. 6. — Ratio of the gas stream to the front disk spectra (see Fig. 4) as a frinction of radius. The notation 
is similar to that of Fig. 5. 
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extended region (e.g., Home et al. 1994). 

In order to investigate the emission along the 
gas stream trajectory, we calculate the ratio of the 
spectrum of the gas stream to the front disk sym- 
metric component at the same radius (Fig. 6). The 
result reveals that the spectrum of the gas stream 
is typically 1 order of magnitude brighter than 
the disk spectrum in the intermediate and mainly 
in the outer disk regions (0.3 < R/Rli < 0.7) 
and this difference is significant at the 3 a level. 
This result suggests the occurrence of gas stream 
overflow. This is in agreement with the results 
of Baptista & Bortoletto (2004), who found clear 
evidence of gas stream overflow in B band eclipse 
maps of V2051 Oph in quiescence. 

In contrast to what is observed in the disk spec- 
tra, in which the lines appear in emission at all 
radii, in the ratio of the gas stream to the disk 
spectrum the CIV, Si III A1892 and Balmer lines 
appear in absorption. This suggests the presence 
of matter outside of the orbital plane (e.g., chro- 
mosphere -|- wind) or that the disk has a non- 
negligible opening angle (flared disk). Horne et 
al. (1994) observed strong CIV emission out of 
the orbital plane in the dwarf nova OY Car in 
quiescence and suggested that this emission may 
arise from magnetic activity on the surface of the 
quiescent disk (Horne & Saar 1991) or from reso- 
nant scattering in a vertically-extended region well 
above the disk. 

The slope of the continuum, the strength of the 
Balmer jump and the line intensities give a wealth 
of information about the physical conditions in the 
emitting gas. It is clear that simple blackbody and 
even stellar atmosphere models arc not adequate 
to describe the strong line emission spectra of the 
V2051 Oph accretion disk. A more quantitative 
analysis of the spatially-resolved spectra demands 
optically thin disk models. This is beyond the 
scope of the present work and will be the subject 
of a forthcoming paper (A. Zabot et al. 2006, in 
preparation). 

3.3. The emission lines 

In this section we analyze the radial behavior 
of selected emission lines. Fig. 7 shows the radial 
intensity distributions for the lines and the adja- 
cent contimmm (top), the radial distributions of 
the net line emission (second from top), the ra- 



dial run of the equivalent width (EW; third from 
top), and the radial run of the full width at half- 
maximum (FWHM; bottom) for the H^, H7, He I 
4471 and He II 4686 lines. The diagrams were com- 
puted for the disk "back" region to avoid contam- 
ination by gas stream emission and, as in Section 
3.2, we used the symmetric disk-emission compo- 
nent in this analysis. We remind the reader that 
our line eclipse maps comprise the line emission 
plus any subjacent continuum contribution. The 
line distributions were obtained from the average 
of all eclipse maps along the line, while the contin- 
uum distributions were obtained from the average 
of the nearest continuum maps on both sides of 
each line. The net line emission distributions were 
computed by the subtraction of the adjacent con- 
tinuum from the corresponding line distributions. 

As expected, the intensity in the line maps is 
stronger than in the adjacent continuum (indicat- 
ing that the lines appear in emission), except in 
the innermost disk regions, where the intensity of 
the continuum reaches that of the lines. For H(5, 
He I 4471 and Hell 4686, the intensity in the line 
is equal or smaller than that of the adjacent con- 
tinuum in the innermost annulus (indicating that 
the lines are in absorption in the white dwarf at 
the disk center) and were not plotted in the cor- 
responding net line emission panels. The slope of 
the continuum distributions is steeper than that 
of the lines. The net line distributions decrease in 
strength with increasing disk radius. The lines are 
in emission at all disk radii with a radial depen- 
dency I oc i?-i-78±o.06 rpj^jg jg steeper than the 
empirical / cx R~^-^ law derived by Marsh et al. 
(1990) assuming a Keplcrian distribution of veloc- 
ities for the line emitting gas in the dwarf nova U 
Geminorum. 

As a consequence of the radial behavior of the 
line and contimmm intensity distributions, the 
lines are relatively weaker in the inner disk re- 
gions and their EW increases with increasing disk 
radius. H(5 and H7 arc the dominant lines with 
an EW ~ 300 A at the disk edge (=i 0.5i?Li), 
while He I 4471 and Hell 4686 have EW ~ 100 A 
at the same radius. In the innermost disk region 
(i? < 0.1 i?Li) the EW becomes negligible (or neg- 
ative) because the continuum intensity reaches (or 
exceeds) the line intensities. Because of this, the 
FWHM in this region is also negligible and is not 
plotted in the corresponding FWHM panels. 
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Fig. 7. — Top: Radial intensity distributions of selected optical lines (filled symbols) and of the corresponding 
adjacent continuum (open symbols). Dashed lines indicate the la limit in each case. Dotted lines mark 
the outer disk radius of i? = 0.56i?Li (Baptista et al. 1998). Second from the top: Net line emission radial 
distributions. Dotted lines show the radial dependency for a slope / oc R~^ '^^. Third from top: EW as a 
function of radius. Bottom: FWHM of the lines as a function of radius. Dashed lines show the law v cx 
expected for a gas rotating in Keplerian orbits. 
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H(5 and He II 4686 show FWHM values compa- 
rable to those expected for a gas rotating in Keple- 
rian orbits around an Mi = 0.78 M© white dwarf, 
but the slope of the radial distribution is flatter 
than the v oc R~^/^ law. Because of the large un- 
certainties in the FWHM values, the distributions 
arc still consistent with the Keplerian expectation. 
On the other hand, the H7 and He I 4471 distri- 
butions are clearly different from the Keplerian 
expectation. Both lines show sub-Kcplcrian ve- 
locities (at the 2 — 3 a confidence level) . While the 
slope of the H7 distribution is consistent with the 
V oc R~^^'^ law, He I 4471 shows a flat distribution, 
with velocities of 1/2 FWHM ~ 600 - 700 km s"! 
at all radii. 

3.4. The uneclipsed component 

The uneclipsed component was introduced in 
the eclipse mapping method to account for the 
fraction of the total light that is not coming from 
the accretion disk plane (e.g. light from the sec- 
ondary star or from a vertically-extended disk 
wind; Rutten et al. 1992). 

Fig. 8 shows the spectrum of the uneclipsed 
component, as well as its fractional contribution, 
function of wavelength. We estimated the 
fractional contribution of the uneclipsed compo- 
nent to the total flux by computing the ratio of 
the uneclipsed light to the average out of eclipse 
level at the corresponding passband. 

The uneclipsed component is dominated by the 
Balmer continuum emission, strong Balmer and 
He emission lines in the optical, and resonant emis- 
sion lines in the UV, suggesting that this light 
arises in a hot, vertically extended, optically thin 
gas region. Lya is the strongest line in the un- 
eclipsed spectrum (~ 30 per cent of the total 
light at that wavelength). Nevertheless, most of 
this contribution is of geocoronal origin and not 
intrinsic to V2051 Oph. The remaining strong 
lines in the optical and UV spectra, as well as the 
Balmer continuum, account for ~ 10 per cent of 
the total light at the respective wavelength. The 
other continuum regions in the uneclipsed spec- 
trum yield negligible fractional contributions to 
the total light. This indicates that the contin- 
uum emission comes mostly from the orbital plane 
while part (~ 10 per cent) of the line emission 
arises in a vertically extended region above the 
disk. 



4. Discussion 

4.1. Evidences for chromospheric emission 

Our continuum eclipse maps, particularly in the 
UV (Fig. 2), show a front-back brightness asym- 
metry. This becomes clear when we compare the 
spectrum of the front and back disk sides (Fig. 5): 
the front side of the disk is brighter than the back 
side. A similar result was found by Vriclmann et 
al. (2002). They discussed a few possibilities (non- 
negligible disk opening angle, warped disk, and 
bulge caused by impact of gas stream) but found 
no compelling explanation for the effect. 

Our front and back disk spectra were computed 
from a symmetric disk component (see Section 
3.2), which traces the baseline of the radial in- 
tensity distribution and removes the contribution 
of asymmetric brightness sources. Furthermore, 
the spectrum of the gas stream region is markedly 
different from the disk spectrum at intermediate 
and large radii {R > 0.3 Rli)- We are, therefore, 
confident that the observed front-back asymmetry 
is not an artifact caused by azimuthal smearing of 
bright spot emission, but it is a real effect. 

Vrielmann et al. (2002) discussed the possibil- 
ity that the observed asymmetry could be caused 
by a non-negligible opening angle (a) in the ac- 
cretion disk, but they discarded this hypothesis 
because it would lead to an enhancement of the 
emission from the back side (seen at the lower ef- 
fective inclination ig// = i — a) in comparison to 
the front disk side (seen at a higher effective in- 
clination if.ff = i + a). This argument is only 
correct if the disk atmosphere shows a stellar-type 
temperature decreasing with vertical height. If the 
disk has a chromosphere, in which the temper- 
ature actually increases with vertical height, the 
effect will be the opposite. Because the front disk 
side is seen at higher angles, the emerging photons 
will sample the uppermost (and hotter) chromo- 
spheric layers and the resulting spectrum will be 
brighter. The back disk side is seen through a 
lower inclination, allowing the emerging photons 
to come from deeper (and cooler) layers and lead- 
ing to fainter intensities. Preliminary fits of the 
spatially resolved spectrum with LTE atmosphere 
models lead to effective temperatures for the front 
disk spectrum that are systematically higher than 
the temperatures of the back disk spectrum at the 
same radius (A. Zabot et al. 2006, in preparation). 
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Fig. 8. — Spectrum of the uneclipsed component. Top: total contribution. Bottom: fractional contribution. 
The notation is similar to that of Fig. 5. 



confirming the above scenario. 

Therefore, the front-back asymmetry observed 
in the accretion disk of V2051 Oph can be inter- 
preted as evidence of chromospheric emission from 
a disk with a non-negligible (but probably small) 
opening angle. This chromosphere is the site of 
the emission lines and is probably also responsible 
for the veiling of the continuum emission from the 
compact underlying white dwarf and bright spot. 
It may also account for the Fell absorption fea- 
tures that become more pronounced with increas- 
ing radius in the disk back side, because the pho- 
tons arising from these regions travel an increas- 
ingly larger path across the vertically-extended 
disk chromosphere before leaving the binary and 
reaching the observer. 

4.2. The white dwarf spectrum and the 
distance 

Here we fit white dwarf atmosphere models to 
the white dwarf spectrum in order to estimate the 
temperature of the primary star and the distance 
to the object. We adopted a white dwarf radius 
of RwD = 0.0103 i?0 and a Roche lobe size of 



Rli = 0.422 i?0 (Baptistaet al. 1998). Therefore, 
our eclipse maps (see Section 2.3) have a scale fac- 
tor of 0.0392 Rli pixel~^ . Since the primary com- 
prises a diameter of 0.0488 the central pixel 
of the eclipse map is fully contained in the white 
dwarf surface. We thus obtain a good estimate of 
the white dwarf spectrum by extracting the flux of 
the central pixel of the eclipse maps at each wave- 
length. These fluxes are multiplied by a factor 



A WD _ / 0.0244Rli \ tt 
Apix ~ V0-0392RLiy' cos^ 



(2) 



to scale the spectrum to the effective area of the 
white dwarf. The resulting spectrum shows a con- 
tinuum filled with broad and shallow absorption 
lines plus a Balmer jump in absorption, resem- 
bling that of a DA white dwarf (Fig. 9). 

We employed a grid of DA white dwarf spectra 
with 6500 K < Two < 20000 JsT and log g = 8 
(D. Koestcr 2000, private communication) in or- 
der to fit the V2051 Oph white dwarf spectrum. 
Since the DA white dwarf models only account for 
the continuum level and the HI lines, we masked 
the other line regions for the fitting procedure. 
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This includes removing the spectral region above 
2300 A in the UV spectrum and below 3400 A in 
the optical spectrum, to avoid contamination by 
absorption bands due to Fell. We considered two 
possibilities: (i) the inner disk is opaque and the 
visible part of the white dwarf surface is the pro- 
jected area of the upper half hemisphere above the 
opaque disk, (ii) the white dwarf surface is fully 
visible through an optically thin disk. 

Our best fit leads to a temperature of Twd = 
9500 1 1900 for the white dwarf with a distance 
of d = 67 Pc if the inner disk region is opaque 
and = 92 +3° pc if the inner disk is optically thin. 
Both the fitted temperature and the resulting dis- 
tance are significantly different from previously re- 
ported values. 

Catalan et al. (1998) made a preliminary fit 
to the G160L extracted white dwarf spectrum 
(assuming a white dwarf atmosphere model plus 
an intervening cool gas layer to account for the 
Fe II absorption bands) to find a white dwarf tem- 
perature of 15000 /C. Their solution provides a 
reasonable fit in the UV but underestimates the 
white dwarf flux in the optical by a factor of 5. 
The white dwarf must be cooler (and, therefore, 
the distance must be smaller) to match the slope 
of the combined UV-optical spectrum. On the 
other hand, Vrielmann et al. (2002) fitted white 
dwarf atmosphere models to UBVRI photometric 
measurements to find Twd = 19600 (if only 
the upper half of the white dwarf is visible) or 
TwD = 15000 if (if the white dwarf is fully visi- 
ble) and d = 146 pc. However, their result largely 
overestimates the white dwarf contribution in the 
UV. For comparison, Fig. 9 shows DA white dwarf 
models for the temperatures and respective dis- 
tances obtained by the works mentioned above. 

As a first step in trying to miderstand the dif- 
ference, we investigate the possibility that some 
systematic error in the eclipse-mapping analysis 
(e.g. wavelength dependent blurring of bright cen- 
tral sources) is affecting the slope of the extracted 
white dwarf spectrum. To test this hypothesis, 
we extract the white dwarf spectrum by measur- 
ing the height of the jump in flux at white dwarf 
egress in the set of narrow band UV and optical 
light curves used for eclipse mapping. This pro- 
vides an estimate of the white dwarf flux as seen 
during egress and is quite similar to the proce- 
dure used by Catalan et al. (1998) to extract the 



UV white dwarf spectrum. The resulting spec- 
trum is consistent with that of Fig. 9 under the 
uncertainties. We are therefore confident that the 
red slope of the white dwarf spectrum of V2051 
Oph at the time of the HST observations is real. 
Thus, the main point is that, it is not possible to 
model the observed white dwarf spectrum with a 
hot {TwD > 14000 /sT) DA atmosphere model. 

The origin of the discrepancy between our re- 
sult and the previous analyses may be in the; re- 
stricted spectral coverage of the data used in those 
works. Catalan et al. (1998) used only the UV 
part of the spectrum (1150 — 2500^1) to fit the 
white dwarf temperature. It is possible to match 
the observed UV intensities with a DA spectrum 
of 15000 K for a distance of d = 184 pc (Vrielmann 
et al. 2002) at the expense of a large mismatch in 
the optical. 

The slope of a DA white dwarf spectrum is 

very similar in the optical (3500 — 5000 A) for a 
large range of temperatures. Therefore, the opti- 
cal spectrum alone does not provide a strong con- 
straint to the white dwarf temperature. We tested 
this by fitting white dwarf models only to the op- 
tical part of the spectrum. This exercise shows 
that it is possible to fit temperatures in the range 
6500 - 16500 K with a only 20% larger than 
the of the best-fit model in this case, T\yd = 
11500 if. This lack of sensitivity affects the white 
dwarf temperature determination of Vrielmann et 
al. (2002), since they used only optical photomet- 
ric data in their analysis. Moreover, they mistak- 
enly took the value Rwd — 0.0244 Rq as the white 
dwarf radius, while the value quoted by Baptista 
et al. (1998) is Rwd = 0.0103 i?,Q. To compen- 
sate for the larger radius assumed for the white 
dwarf, their fit had to increase the distance by a 
similar amount (a factor of ^ 2.3). 

Steeghs et al. (2001) used optical spectroscopy 
(4000-6800 A) to find that a (15000 K) blackbody 
was a good match to the slope of their extracted 
white dwarf spectrum. Their result faces the same 
problem discussed above; the slope of the optical 
continuum is very similar for a large range of tem- 
peratures. 

It is important to notice that at the time of 
our observations V2051 Oph was in an unusual 
low brightness state, in which mass accretion may 
have been considerably reduced (Baptista et al. 
1998). With reduced (or absent) accretional heat- 
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wavelength (A) 

Fig. 9. — Spectrum of the V2051 Oph white dwarf (histogram) and DA white dwarf model fits. The dot- 
dashed Une corresponds to a model with T — \%QQQK^d — 146 pc (Vrielmann et al. 2002), and the dashed 
Hne correspond to a model with T — lhOOQK,d — 184pc (Catalan et al. 1998). The best fit model is 
displayed as a solid line. The bottom panel shows the S/N of the the white dwarf spectrum. 
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ing, it is possible for the white dwarf to cool down. 
We therefore expect that our fitted white dwarf 
temperature is lower than that during a normal, 
quiescent state. 

Our larger spectral coverage, combining optical 
and UV data, allows a better determination of the 
slope of the white dwarf spectrum, and, therefore, 
its temperature and distance. Our results show 
that the distance to V2051 Oph is smaller than 
previously found. 

5. Conclusions 

The main results of our spectroscopic study of 
V2051 Oph in a faint brightness state during 1996 
can be summarized as follows: 

1. The presence of white dwarf and bright-spot 

strong emission in the continuum maps and 
their absence in the line maps, coupled with 
the significant extra absorption in the spec- 
tra of the back disk side, arc indications 
of the vertical extension and large optical 
depth of the gas from which the lines origi- 
nate. 

2. Distinct emission along the stream trajec- 
tory suggests the occurrence of gas stream 
overflow. 



7. The front disk spectrum is systematically 
brighter than the back disk spectrum at the 
same radius. This can be explained in terms 
of chromosphcric emission (higher tempera- 
tures at the uppermost atmospheric layers) 
from an accretion disk with a non-negligible 
opening angle (limb brightening effect). 

8. We fit stellar atmosphere models to the ex- 
tracted white dwarf spectrum to find a tem- 
perature TwD = 9500 +11'^'^ K and a dis- 
tance of d = QT^'^pc (if the inner disk is 
opaque) or d = 92 pc (if the inner disk is 
optically thin). 
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3. Spatially resolved spectra show that the lines 
are in emission at all disk radii. The Balmer 
decrement becomes shallower with increas- 
ing radius. 

4. The FWHM of the emission lines differs from 
that expected for a gas in Keplerian rotation 
and the line intensities drop with a radial 

dependency of / oc i?-i-78±o.06^ 

5. The spectrum of the uneclipsed light is dom- 
inated by strong emission lines and a Balmer 
jump in emission, indicating origin in a hot, 
vertically extended, optically thin gas region 
above the disk. The strongest uneclipsed 
lines contribute ~ 10 per cent of the total 
flux. 

6. Strong absorption bands, possibly due to Fe 
II, are seen in the spectra of the back disk 
side, suggesting it arises from absorption by 
an extended gas region above the disk. 
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A. Reconstruction of maps from light curves with incomplete pheise coverage 

Here wc! address the reliability of eclipse-mapping brightness reconstructions for the case of light curves 
with incomplete phase coverage, such as those of run HI. 

In order to simulate the structures observed in the optical maps, we created artificial eclipse maps of 
51 X 51 pixels with two compact Gaussian spots, one at the disk center (the white dwarf) and another near 
the impact region of the gas stream with the disk (the bright spot), on top of a low- intensity, broad Gaussian 
background to simulate a faint accretion disk. The artificial map is shown in Fig. 10. 

We simulated the eclipse of the artificial map using the geometry of V2051 Oph {i = 83° and q = 0.19, 
Baptista et al. 1998) to create light curves with added Gaussian noise reproducing the S/N of the V2051 
Oph data {S/N = 10). Since the optical light curves of V2051 Oph have incomplete phase coverage, we 
truncated the synthetic light curves to reproduce the set of orbital phases of the optical run. Fig. 10 shows 
light curves for a complete phase coverage and for the set of phases of the V2051 Oph optical data. We 
applied eclipse mapping techniques to the two light curves and compared the results with the original map. 

The ability of the eclipse-mapping method to construct a two-dimensional surface brightness map from 
a one-dimensional eclipse light curve stems from the fact that ingress and egress Hrdic;s intersect with each 
other at relatively large angles (Horne 1983). Fig. 11 illustrates the eclipse geometry of V2051 Oph for the 
incomplete set of phases of run HI. Arches drawn on the eclipse map connect points of constant ingress 
and egress phases. The pair of thick lines shows the ingress/egress arches for the first phase of the light 
curve {(j) = -1-0.0194). Each pair of arches outlines a region of the disk that is being occulted by the shadow 
of the secondary star at a given phase (for visualization purposes, Fig. 11 shows only one of every three 
phases/arches in the incomplete light curve). 

The effect of incomplete phase coverage is to reduce the information available to reproduce structures 
in parts of the accretion disk. For the light curve with incomplete phases, about 40 % of the disk surface 
is not covered by ingress or egress arches (Fig. 11). There is no information in the eclipse shape about 
the brightness distribution of the disk regions not sampled by the grid of arches. Therefore, the disk side 
approaching the secondary star (the lower hemisphere of the map in Fig. 10) will not affect the shape of the 
eclipse light curve with the incomplete set of phases and its brightness distribution will not be recovered by 
the eclipse-mapping method. 

Fig. 10 shows that the structures are reasonably well reconstructed with the light curve of full eclipse 
phase coverage. As a consequence of the intrinsic azimuthal smearing effect of the eclipse mapping technique 
(see Baptista 2001 for more details) the asymmetric spot appears smeared into a crescent-shaped structure. 
Because the light curve with incomplete phase coverage lacks most of the ingress phases, the technique is not 
able to properly reconstruct structures located in the disk region that approaches the secondary. However, 
for this V2051 Oph simulation the main brightness structures are in the disk center and in the bright spot 
region (the disk hemisphere moving away from the secondary star), and the reconstruction is reasonably 
good and quite similar to that obtained with the complete eclipse phase coverage. The radial intensity 
distributions of the eclipse maps from the complete and incomplete phase coverage light curves (not shown) 
are very similar, with only small differences that are within the uncertainties of the distributions. 

In summary, the simulations show that, although the eclipse mapping method is not able to recover 
structures in disk regions for which there is no information in the eclipse shape, brightness distributions with 
structures in the disk center and in the bright spot region can still be detected in eclipse maps derived from 
light curves with an incomplete phase coverage such as that of run HI, in spite of their relatively low S/N. 
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Fig. 10. — Left: Synthetic light curves with S/N=10 for the full phase coverage (top) and for the phase set 
of the optical run (bottom). Right: Original map and reconstructions obtained from the light curves in the 
left panels. The notation is the same as in Fig. 2. 
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Fig. 11. — Eclipse geometry of V2051 Oph for the set of phases of the optical run. Arches (solid lines) 
connect points of constant ingress and egress phases. To avoid overcrowding the figure, only one phase/arch 
of every three is shown. The pair of thick lines shows the ingress and egress arches for the first phase of the 
light curve (0 = +0.0194). Dotted hues depict the primary Roche lobe, a disk of radius 0.56 Rli and the 
ballistic stream trajectory. 
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